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ABSTRACT 

While all the cosmological observations are carried out on a light-cone, the null 
hypersurface of an observer at z = 0, the clustering statistics has been properly defined 
only on the constant-time hypersurface. We develop a theoretical formulation for a 
two-point correlation function on the light-cone, and derive a practical approximate 
expression relevant to the discussion of clustering of high-redshift objects at large 
separations. As a specific example, we present predictions of the two-point correlation 
function for the Durham/AAT, SDSS and 2dF quasar catalogues. We also briefly 
discuss the effects of adopted luminosity function, cosmological parameters and bias 
models on the correlation function on the light-cone. 
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1. Introduction 

Improving magnitude-limits of astronomical surveys naturally increases the fraction, and 
therefore the weight, of selected objects towards higher redshifts in the entire sample. In fact 
discussion of clustering of objects at z = (1 ~ 3) is becoming fairly common, including the 
Lyman-break galaxies (Steidel et al. 1998; Jing & Suto 1998), X-ray selected AGNs (Carrera et al. 
1998), the FIRST survey (Magliocchetti et al. 1998), and up-coming 2dF (2 degree Field; Boyle 
et al. 1998) and SDSS (Sloan Digital Sky Survey) QSO surveys. The clustering statistics of such 
high-z objects provides several important pieces of cosmological information albeit in a rather 
complicated manner such as linear and nonlinear evolution of mass density fluctuations (Hamilton 
et al. 1991; Jain, Mo, & White 1995; Peacock & Dodds 1994,1996), evolution of object-dependent 
bias (Mo & White 1996; Jing 1998; Fang & Jing 1998), and redshift-space distortion (Kaiser 1987; 
Hamilton 1997; Ballinger, Peacock & Heavens 1996; Matsubara & Suto 1996). 

Another important but less often discussed point is the light-cone effect, that is, such 
cosmological observations are feasible only on the light-cone hypersurface defined by an observer 
at z = 0. In the case of an angular two-point correlation function of high-z objects, the light-cone 
effect can be included relatively easily (e.g., Lahav et al. 1997; Yamamoto &; Sugiyama 1998), 
because all the positions of objects are projected on a two-dimensional sphere on the sky. In 
discussing the spatial two-point correlation function of these objects, however, the light-cone effect 
hampers any attempt to distinguish the scale-dependence of clustering in the survey volume from 
the intrinsic redshift evolution (e.g., change of the mean number density of the objects considered). 

Some aspects of the light-cone effect have been already discussed by Matarrese et al. (1997), 
Matsubara, Suto & Szapudi (1997), and Nakamura, Matsubara & Suto (1998). In those papers 
the light-cone effect is taken into account by integrating over the line-of sight convolved with the 
selection function. Matarrese et al. (1997), for instance, adopted an expression: 
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where dN/dz is the observed differential redshift number count of the objects, b e g(z) is the effective 
biasing factor for the objects, and £(r, z) is the mass two-point correlation function at a comoving 
separation r. 

While the above expression looks physically reasonable, it is important to properly define 
the two-point correlation function on the light-cone hypersurface and then to derive useful and 
practical expressions with clarifying the underlying approximations and assumptions. Even apart 
from such a theoretical motivation, equation (|I|) involves a double integration over the redshift 
which is not easy to evaluate numerically. Thus it is desirable to find a more practical approximate 
formula, if any. 

The primary purpose of this paper is to present a rigorous theoretical formulation to define 
and compute a two-point correlation function of cosmological objects which explicitly takes 
into account the light-cone effect for the first time (§2). We propose an expression (eq.^J) as 
a practically useful formula relevant to the discussion of clustering of high-z objects at large 
separations (§3). As a specific example of high-redshift objects on the light-cone, we first consider 
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the Durham/AAT QSO samples (Shanks & Boyle 1994; Croom & Shanks 1996) in §4.1. We 
compute the corresponding two-point correlation function on the basis of our formula (eq.p^j), 
which turns out to be in good agreement with the results of Matarrese et al. (1997) which 
adopts equation ([l]). Unfortunately we are not able to make further quantitative comparison with 
their predictions due to the fact that equation (|l|) is hard to evaluate numerically. Therefore we 
illustrate the importance of the light-cone effect by comparing with predictions based on another 
intuitive expression (eq. |24|]) which we propose as a counter example. For that purpose, we start 
from the QSO luminosity function by Wallington &: Narayan (1993), and predict the two-point 
correlation functions corresponding to the future SDSS and 2dF QSO samples (§4.2). Also we 
briefly discuss the dependence on the underlying cosmological parameters (§4.3). 

We should note here that the present paper is rather theoretical at this stage. While 
we propose equation ( f23|) as a well-defined and practically useful expression for the light-cone 
correlation function, the crucial difference shows up only in the regime where the correlation is 
very weak which is hardly probed accurately with the current data like Durham/AAT sample, 
for instance. Furthermore we have not yet successfully included several important effects 
observationally; redshift-distortion and evolution of bias among others. We plan to report the 
further progress in those issues elsewhere in due course (Nishioka & Yamamoto 1999; Yamamoto 
& Suto 1999). Such limitations are discussed in §5 along with our main conclusions of the present 
paper. Throughout this paper we use the units in which the light velocity c is unity. 



2. Two-point clustering statistics on a light-cone hypersurface 

In what follows, we focus on the spatially- flat Friedmann - Robertson - Walker space-time 
for simplicity, whose line element is expressed in terms of the conformal time r/ as 

ds 2 = a 2 (ri) \-drj 2 + dr 2 + r 2 dn 2 {2) ] . (2) 

Since our fiducial observer is located at the origin of the coordinates (ry = r = 0), an object at 
r and r] on the light-cone hypersurface of the observer satisfies a simple relation of r = t]q — rj. 
While this is mainly why we adopt the spatially-flat model, current observations indeed seem to 
favor such a cosmological model (e.g., Garnavich et al. 1998). 

Denote the comoving number density of observed objects at r\ and x = (r, 7) by n(rj,x), then 
the corresponding number density defined on the light-cone is written as 

n LC (r,7) = n(r/ - r,r,7). (3) 

If we introduce the mean observed number density (comoving) and the density fluctuation at 77, 
tiq{j]) and A(ry,x), on the constant-time hypersurface: 

n(??,x) = n Q {rj) [1 + A(ry,x)] , (4) 

equation (||) is rewritten as 

n LC (r, 7) = n (r/o ~r)[l + A(t? - r, r, 7)] . (5) 
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Note that the mean observed number density no(rj) is different from the mean number density of 
the objects 77(77) at rj by a factor of the selection function <f>(r]) which depends on the luminosity 
function of the objects and thus the magnitude-limit of the survey, for instance: 

710(77) = 71(77)0(77). (6) 

When the observed density field of objects on the light-cone, n LC (r, 7) is given, one may 
compute the following two-point statistics: 

1 f (7f2 f f 
X ( R ) = yLC J -4^ J dnrfdn^J dr 2 rid^ 2 n LC (r 1 ,7 1 )n LC (r 2 ,7 2 )5( 3 )(x 1 - x 2 - R), (7) 

where xi = (ri,ri7i) and x 2 = (r 2 ,r 2 7 2 ) and R = |R|, R = H/R, and V LC is the comoving survey 
volume of the data catalogue: 

V LC = f ^ r 2 dr f c?% = ^(r max 3 - r min 3 ), (8) 

with r max = r(z max ) and r m i n = r(z m i n ) being the boundaries of the survey volume. Although the 
second equality as well as the analysis below assumes that the survey volume extends 47r steradian, 
all the results below can be easily generalized to the case of the finite angular extension. 

Substituting equation (||), the ensemble average of an estimator X{R) is explicitly written as 

(X(R)) =U{R) + W(R), (9) 

where 

1 r dn T 



12 



U ( R ) = vlcJ^tJ drir "J dn * J dr2T2 S dn ^ 

xn (r/o - T-iHfao - 7- 2 )5 (3) (xi - x 2 - R), (10) 

and 

= 7/Cc / J drir i J dn ii J dr2r 2 / dn^n (m ~ n)n (vo ~ r 2 ) 

x(A(t7o - ri,ri,7i)A(77o - 7- 2 ,r 2 ,7 2 )W 3) (xi - x 2 - R). (11) 



Consider first W(R) which contains all the information of the clustering. We show in 
Appendix A that in linear theory W(R) reduces to 

W ( R "> = y^wR I Is dridr2rir2 I! M% ~ r j) D i(Vo ~ rj)\ 

x / dkk 2 P(k)j (kR)b(k; 7/o - n)b(k; T/o - r 2 ), (12) 



where D\{rj) is the linear growth rate normalized to unity at present (770), P(k) is the power 
spectrum of the mass fluctuations at 770, b(k;rj) is the /c-dependent linear bias factor, j'o(x) is the 
spherical Bessel function of the 0-th order, and S denotes the region \r% — r 2 | < R < 7*1 + r 2 . 
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Since we are generally interested in the case of R <C r max , we can use the approximation: 

/■rmax rR+ri 

dr\dr 2 — / dr\ / dr 2 . (13) 

S Jr min J\R-r 1 \ 

Furthermore using that r/o — r 2 ~ r/o — r\ , equation (|l^) reduces to 

W(R) ~ / r 2 dr [n (?7o - r)Di(% - r)] 2 

* •''"min 

X i / k 2 dkP(k)j (kR) [b(k; m - r)] 2 . (14) 

Therefore in linear theory, we can write 

47r r max o o 

VO ~ J") Source , (15) 

where £(i?; 77) source is the conventional two-point correlation defined on the constant hypersurface 
at the source's position: 

resource = 2^2 / k 2 dkP(k)j (kR) [&(£; r,)^ (r,)] 2 . (16) 



Consider next U(R). Repeating the similar calculation from equation (|il| ) to (|l2|), equation 
(|To| ) reduces to 

U ( R ) = TFTc~n 1 / dridr 2 rir 2 n (ri - ri)n (?7o - t* 2 ), (17) 



and the same approximation from equation (|l^) to ( [To]) yields 



47T /" r max 

"(^ " T^lc / r 2 dr[n (r?o-r)] 2 . (18) 

Thus W(i?) is independent of R for i? -C r max , as expected. 

Although the rigorous derivation in Appendix A is carried out entirely in the framework of 
linear theory, the final result ( |i~5| ) would be valid also for non-linear regimes; replacing equation 
( |i~6| ) by its nonlinear counterpart, one may approximate the correlation in equation (11) as 



^A(r/ - ri,ri,7i)A(?7o - r 2 ,r 2 , 72)) ~ £(ri,7* 2 ) ^ £(-R;f?o - Hi), (19) 

where = (7*1 + r 2 )/2, i? = |ri — r 2 |, and ^{R;rj) is the two-point correlation function on 
an equal-time hypersurface. In equation (|l9|) , we have neglected the angular dependence in 
(A (770 — 7*i, 7*1, 71) A (770 — r 2 ,r 2 ,7 2 )) and regarded it as a function of 7*1 and r 2 . If this is the case, 
the similar calculation to derive equation ( |l~7l) yields 



W{R) = yLC~R y y^ dridr 2 rir 2 n (rio - ri) 710(770 - r 2 )£(i?; r/ - r h ). (20) 
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Again, with the same approximation from equation (12) to (|l5|), we have 



W(R) 



Ait 
yLC 



r 2 dr [n (rj - r)] £(R; rj -r), 



(21) 



which is identical to equation (15) except that £,(R;r]o — r) is not restricted to linear theory 
prediction. Unfortunately we have not been able to justify the validity of the approximation ( |l9| ) 
strictly, but we expect that the above qualitative and intuitive argument supports it. Note that 
this approximation is also implicitly assumed in equation (|l]). More importantly the resulting 
non- linear effect, which is taken into account in this way in equation (|2~I|), is small as long as the 
QSO clustering at R £ l/i -1 Mpc is considered as we will explicitly show in §4. 



3. Definitions of the two-point correlation function on the light-cone 

In the previous section, we derived an approximate expression for the two-point statistics 
(X(R)) and W(i?) on the light-cone. The remaining task is to define the corresponding two-point 
correlation function on the light-cone. A straightforward application of equation (||) implies the 
definition: 

LC _ {X(R))-U(R) _ W(R) 
U {R) = U{R) " U{R) ■ (22) 

Substituting equations (|l^) and (18), the above definition reduces to 

/ rmax 2 2 

drr n (r/ - r) £(R; rj - r) Source 

Ci ( "> - • (23) 

drr 2 no(r]o — r) 2 

^min 

Another possibility, which might look more similar to a conventional pair-count estimator adopted 
in analyzing an observational map on a constant-time hypersurface (z ~ 0), is 



& {R) = 7^ = 7^ ' ( 24 ) 



where n LC is the mean number density on the light-cone hypersurface and (n LC ) denotes its 
ensemble average: 

n LC = — r / r 2 dr f dfy n LC (r,j), (25) 



(" LC ) = TTTc I r 2 drn (vo-r). (26) 



yLC 
yLC 

Two-point correlation functions on the light-cone can be computed from the given catalogue 
of objects according to either theoretical definition in a fairly straightforward manner, although 
one has to assume a set of cosmological parameters a priori to translate the observable coordinates 
to the comoving ones; first average over the angular distribution and estimate the differential 
redshift number count dN/dz of the objects. Second distribute random particles over the whole 
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sample volume so that they obey the same dN/dz. Then the conventional pair-count between 
the objects and random particles yields X(R) (although not (X(R)}, of course), while U(R) is 
estimated from the pair-count of the random particles themselves. Finally n LC can be computed 
as J 2 Zmax dz{dN/dz) with an appropriate normalization constant. In addition, when the observable 
comoving number density no(r]) is independent of time, ^J^(R) and ^^{R) become identical and 
reduce to 

/f max 
drr £(R;rfy - r) SoU rce 
ha-w = srw = — • (27) 



max 

drr 



If the correlation function of objects does not evolve, i.e., £(R;r]o — ^Source = £(R] ?7o)source> 
however, equation (|^) readily yields 

Za°(R) = £(#! ^Source (28) 

while equation ( |24|) does not. In fact, as is shown in § 4.2, ^ c (i?) does not vanish even where 
the real spatial correlation is negligible. This is an example of a non-physical effect due to the 
contamination by the intrinsic evolution of mean number density of objects on the light-cone 
hypersurface. 

From a theoretical point of view, one might expect that a more accurate expression on the 
light-cone should involve the double integrals with respect to z (e.g., eq. Q) since the spatial 
correlation along the line-of-sight should manifest itself in the correlation between different 
redshifts. While this argument is correct in a strict sense, we have shown that our expressions 
without the double integral already take into account the leading term properly. Thus unless 
one needs a next-order correction, our approximate formulae are practically sufficient, and better 
suited for numerical evaluations, compared with equation ([l]). 



In summary, we propose equation (p2|), rather than equation (|24|), as the best estimator of the 
two-point correlation function on the light-cone hypersurface which captures the true physics of 
correlation due to clustering. Nevertheless, it is instructive to compare with the latter and discuss 
their difference in some detail in order to understand some aspects of the light-cone effect, as we 
show below. 



4. Application: two-point correlation functions of high-redshift quasar samples on 

the light-cone 

4.1. Durham/AAT QSO sample 

First we examine in details the extent to which equations (|2^), and (|24] ) lead to 
quantitatively different predictions. For this purpose, we compute the correlation functions of 
Durham/AAT QSO sample (Shanks & Boyle 1994; Croom & Shanks 1996) in the Einstein - de 
Sitter universe, following Matarrese et al.(1997). In this case, a{r\) = (r]/r)o) 2 , and the present 
conformal time t/q is 2/Hq where Hq = lOO/i km/s/Mpc is the Hubble constant. 

According to Matarrese et al. (1997), we adopt a simple assumption that the quasar 
correlation function is given by the mass two-point correlation function at that epoch multiplied 
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by the bias factor which depends on the mass of the hosting dark matter halo (e.g., Mo & White 
1996; Mo, Jing, & White 1997; Jing 1998; see also Fang & Jing 1998). In this case, Matarrese 
et al. (1997) found that the effective bias b e g(z) integrated over the halo mass larger than M m ; n is 
well-fitted to 

b eS (z) = 0.41 + [b e{i (z = 0) - 0.41] (1 + zf, (29) 

and obtained the values of b e s(z = 0) and (3 for given M m j n in the standard cold dark matter 
(CDM) model (see their Table 1). As for the mass two-point correlation function, we use 
the nonlinear fitting formula by Peacock &: Dodds (1994,1996) for the CDM power spectrum 
normalized to the cluster number counts (Kitayama & Suto 1997; Kitayama, Sasaki & Suto 1998). 

We adopt the polynomial fit of Shanks & Boyle (1994) to the differential number count of the 
Durham/AAT QSO sample (the corresponding B-band limiting magnitudes in several fields range 
20.12 to 21.27): 

(^p) = 2.738 - 57.573z + 341.720z 2 - 548.737z 3 + 408.908z 4 - 145.698z 5 + 19.989z 6 (30) 

\ dZ J Durham 

for 0.3 < z < 2.2 (392 QSOs in total) which is plotted in Figure |l|a. We translate the count N(z) 
to the observed mean number density no(r) according to 

dN 1 dz 

»W) oc (31) 

dz r z {z) dr 

where the proportional constant is not necessary in our formulae. In this case we found that both 
equations (23) and ( p4| ) yield almost indistinguishable results. Our results are plotted in Figure |2| 



along with those of Matarrese et al. (1997). It is clear that the three expressions for the light-cone 
correlation functions result in negligible difference, at least comparing with both the statistical 
errors from the Durham/AAT QSO sample and the theoretical uncertainty of the evolution of 



bias. In the next subsection, however, we show that equations (|23| ) and (24) exhibit significantly 
behavior at large R where the correlation is weak. Unfortunately we are not able to make further 
comparison with equation (|l[) proposed by Matarrese et al. (1997) since we are not sure how to 
deal with the double integration precisely. 



4.2. 2dF and SDSS QSO samples 

It seems premature to draw any cosmological conclusions from the comparison with the 
currently available data samples. Therefore we apply our formulae, ^\ C (R) and £g C (ii), to predict 
the correlation functions for the ongoing SDSS and 2dF QSO catalogues. In this case we need 
a redshift dependent QSO luminosity function which is not well-established, especially at higher 
redshifts. In what follows we adopt the B-band quasar luminosity function according to Wallington 
& Narayan (1993; see also Boyle, Shanks & Peterson 1988 and Nakamura & Suto 1997). To be 
specific, for 0.3 < z < 3 

*(M B) Z) = 10 o.4(a+l)[M B -M*( 2 )] + 10 0.4(/?+l)[A/ B -Af*( Z )] ' ^ 

M£(z) = M^-2.5fc L log(l + z), (33) 
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with = -20.91 + 51og/i, k L = 3.15, a = -3.79, (3 = -1.44, = 6.4 x 10^ 6 /i 3 Mpc 3 . For 
z > 3, they adopt 

*(M B) *) = 10 o. 4 (a+l)[M B -M*( 2 )] + 10 0.4(/3+l)[A/ B -M*( 2 )] ' ( 34 ) 

M£(z) = Mg — 2.5&L log 4 + (35) 

with A = (z — 3) log 3.2, B = 2.5^4/(a — (3). To compute the B-band apparent magnitude from a 
quasar of absolute magnitude Mb at z (with the luminosity distance d^(z) = (2 / Hq)(1 + z — y/1 + z) 
in the Einstein - de Sitter universe), we applied the K-correction: 

B = M B + 5 log(d L («)/10pc) - 2.5(1 - p) log(l + z) (36) 

for the quasar energy spectrum L v oc v~~ p (we use p = 0.5). We adopt the B-band limiting 
magnitudes for 2dF and SDSS QSO samples as B\ im = 20.85 and 20, respectively. The redshift 
distribution of QSOs predicted from the above luminosity function is plotted in Figure |]b with 
the B-band limiting magnitudes B\ im = (18 ~ 21). 

We show our £^(R) and ^ C (R) for SDSS and 2dF QSO samples together with the nonlinear 
mass two-point correlation functions in Figure ||. Here we adopt the standard cold dark matter 
(SCDM) model, in which 0,q = 1, A = 0.0, h = 0.5, and the amplitude of the CDM density power 
spectrum is normalized as as = 0.56 according to the cluster abundances (e.g., Kitayama & Suto 
1996). Three lines correspond to different threshold mass My im = 10 13 , 10 11 and 10 9 M Q (from top 
to bottom) of the dark matter halos which are supposed to host each QSO (Matarrese et al. 1997). 

We plot the results using the linear fluctuation power spectrum in dashed lines while solid 
lines use the non- linear models of Peacock & Dodds (1996). The nonlinearity becomes important 
only on R ^ l/i _1 Mpc implying that one can safely ignore the nonlinear effect on scales which are 
probed by (sparse) QSO samples in general. The linear and nonlinear correlation functions for 
dark matter (defined on the constant-time hypersurface) are plotted in Figure ||d. 

Panels (a) and (b) adopt the B-band limiting magnitude B\\ m = 20.0 corresponding to the 
SDSS QSO catalogue, while panel (c) adopts B\\ m = 20.85 for the 2dF QSO catalogue. In practice, 
those two predictions are very similar if we use £^ c (i?). On the other hand, £^ C (R) behaves 
quite differently; in fact it levels off at large R, and the asymptotic value depends on the range 
of z (Fig.^b). This is qualitatively explained as follows; since W(R) decreases as R and U{R) is 
independent of R at large R, equations (|T^ ) and (§J) imply that ^ C (R) becomes asymptotically 



K— >oo (n ) 



n j ; i .-. .. 



r dr x r dr[no(r]o — r)]" 

s L ( 37 ) 

r 2 drno(??o — r ) 



The right-hand side of the above equation does not vanish unless n^rf) is constant, and becomes 
larger as no(rf) changes more significantly within the survey volume. This is why ^ C (R) becomes 
significantly different from £^ c (i?) when the survey limit exceeds z ~ 3 where the QSO number 
density begins to decrease substantially (see Fig. |l|a and b). Incidentally this feature was not 
clear in Figure |2| because the Durham/AAT QSO sample is limited to z < 2.2 and also because 
the Figure is plotted in linear, not log, scale. 
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This clearly illustrates the fact that the spatial clustering of objects defined on the light-cone 
could be apparently contaminated and mixed up with the intrinsic evolution of mean density of 
objects and even with the shape of observational selection function. In this respect, our proposal 
of £^ c (i?) is a more reliable and robust definition describing the two-point correlation on the 
light-cone. 



4.3. dependence on the cosmological model parameters 

The next important question is the extent to which different cosmological models lead to 
different predictions of two-point correlations. For this study, we need the QSO luminosity 
function and bias models for the arbitrary cosmological models which are not available at this 
point. So we simply adopt the Durham/AAT QSO redshift distribution function (30) and the bias 
model (^), although the latter is derived in the Einstein - de Sitter universe. 

In Figure EJ, we compare the predictions in the SCDM model against those in the spatially-flat 
low-density CDM (LCDM) model in which tt = 0.3, A = 0.7, h = 0.7, and a 8 = 1.0. Figure 
| plots the amplitude of Ca°( r ) at R = 15/i _1 Mpc function of the redshift of the survey 
limit while keeping B\[ m = 20. Since LCDM has more power on larger scales and the fluctuation 
amplitude erg is larger, the correlation is stronger compared with that in SCDM given the same 
bias model. It should be noted that either model predicts that the QSO correlation amplitude 
increases as z becomes larger unlike the mass correlation which always grows from high to low 
redshifts. This qualitative feature is consistent with the finding of La Franca, Andreani & Cristiani 
(1998), and implies that the evolution of bias dominates the growth of high-z objects in addition 
to the growth rate of the mass fluctuations. Therefore this kind of comparison should yield 
profound cosmological implications on the nature of QSOs, although it would be premature to 
draw any decisive conclusion with the current theoretical understanding of the bias (e.g., Fry 1996; 
Mo & White 1996; Mo. Jing, & White 1997; Fang & Jing 1998; Jing 1998) and statistics of the 
observational samples. 



5. Discussion and conclusions 

Although the clustering of objects at high redshifts has been extensively discussed in the 
literature, the previous theoretical analysis was largely based on a more or less qualitative 
treatment of the light-cone effect (Matarrese et al. 1997; Matsubara, Suto &: Szapudi 1997; 
Nakamura, Matsubara & Suto 1998). For on-going wide and deep surveys like SDSS and 2dF QSO 
surveys, the light-cone effect becomes very important either as a contamination of the real signal 
or as a cosmological probe. 

In the present paper, we developed a theoretical formulation which properly takes account of 
the light-cone effect for the first time. Strictly speaking we were able to derive our main result, 
equation (|l5|), only in linear theory, but we expect that the same expression would be a good 
approximation even in nonlinear regime. In any case the correlations of objects at high redshifts 
are described almost entirely in linear theory for R ^ l/i _1 Mpc (Figs.|| and ||), and in practice the 
expression is guaranteed to be valid on the scales of interest. 



We proposed a well-defined expression, equation (22), for two-point correlation functions on 
light-cone, and compared its implications with those of another possibility (24). As long as the 
Durham/AAT QSO sample is considered, both expressions give almost indistinguishable results 
and in fact they are also in good agreement with the previous proposal by Matarrese et al. (1997). 
Although we are not able to make further quantitative comparison with the latter, this implies 
that our proposal without the double integration over the redshift distribution is more practical 
in making theoretical predictions. Then we applied our expressions and computed the two-point 
correlation functions on light-cone for future SDSS and 2dF QSO samples, and showed an example 
that the light-cone effect could mix up the true spatial clustering of objects and the intrinsic 
evolution of mean density of objects if one uses a native definition like equation (|24T). 



In fact there are several issues which remain to be worked out in the present context. The 
QSO luminosity function and the evolution of bias play central roles in confronting observations 
and predictions of the QSO correlation functions. In this paper, we tentatively adopted the 
expression by Wallington & Narayan (1993) even in LCDM models, although it is relevant only in 
the Einstein - de Sitter model. It is highly desirable that the QSO luminosity function in general 
cosmological models is derived from future observations and becomes available for the theoretical 
analysis. Theoretical approaches to determine the time (and scale) dependence of bias are just in 
the beginning (Fry 1996; Mo & White 1996; Jing 1998; Dekel & Lahav 1999). We did not attempt 
to explore a range of possible bias models, but rather adopt a simple fit by Matarrese et al. (1997) 
as a specific example. Most likely QSO correlation functions from future samples are the most 
straightforward tool to test the several bias models in further detail. In the present paper, we 
have neglected the redshift-space distortion either to the peculiar velocity of the objects (Kaiser 
1987; Hamilton 1997; Nishioka & Yamamoto 1999) or to the geometry of the universe (Matsubara 
& Suto 1996; Ballinger, Peacock &: Heavens 1996). Definitely these effects should be crucial in the 
quantitative comparison to high precision, and we plan to incorporate the effect in future work. 
Nevertheless we hope that the current paper presents a convincing case that the light-cone effect 
should be properly taken into account in analysing the future surveys of high-redshifts objects. 
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Appendices 

A. Calculation of W{R) 
We present an explicit derivation of equation (|l~2^) from equation (11): 

W{R) = / / drir i / d ^7i / dr2r 2 / dtt^ 2 n (r) - ri)n (??o - r 2 ) 

x/A(?7o -n,ri,7i)A(r/ - r 2 ,r 2 ,7 2 )W 3) (xi — x 2 — R) (Al) 



in linear theory. 

First let us expand the number density contrast A(rj,r, 7) as 



poo 

A(r/,r,7)=/ dk J2 &kim{n)y k im(r, 7), (A2) 



J,m 

in terms of the normalized harmonics: 

yfcjm(r,7) = X jU (r)y lm (n f ) J (A3) 

where 

X M (r) = J^kj^kr), (A4) 

V 7T 

and Y/ m (f2^) and ji(x) are the spherical harmonics and the spherical Bessel function, respectively. 
Substituting equation (|A2[) into equation JA1|) , one obtains 

= ^Tc / / drir i J dfl "fi J dr2r 2 J d^2 n o(Vo ~ ri)n (r?o - r 2 ) 
x J dh J dk 2 ( A kihmdVo ~ rx)Al 2hm2 (rio - r 2 )) 

h,mi lifmi 

x^^n)^™,^,)^^)^^^) ^ 3 )( Xl - x 2 - R). (A5) 

In Appendix B, we show that the ensemble average of the mode coefficient in the integrand 
of equation ( |A5| ) reduces to the following function: 

(Afc liimi (77b - ri)Afc 2 ; 2m2 (r/ - r 2 )) = Q{ki,k2,r l ,r 2 ,h)5 hh 5 mim2 , (A6) 
where 5y 2 and 5 mim2 are the Kronecker's delta. In addition, we use the relations 

<^( Xl - x 2 - R) = | d 3 k e-ik-Cxx-x.-R)^ (A?) 

and 

e -&* = 4 ^]T £ (-O'iiW^imCnfc)^^), (A8) 

i m=—l 



14 



and then equation ( |A5| ) becomes 

1 r dn 



W(R) = y^c J J drir i / d ^7i / dr2r 2 / d ^72 n o(% - ri)n (r? - r 2 ) 
X J dki ^2 J dk 2 ^2 Q(h,k2,r 1 ,r 2 ,h)Si 1 i 2 5 mim2 

h,mi i 2 ,m 2 
xX fclil (ri)^ imi (^ 1 )X fc ^ 2 (r 2 )y i * m2 (% 2 ) 



X 



^3 f d 3 k4. Yl (-^^M^W^lW^ 

" ! 7 LiAfi 

x4vr ^ (0 X2 ii 2 (fcr 2 )y L * 2 M 2 (^k)^ 2 M 2 (n f2 



(27T) 



L2M2 

x4vr WSiJ^W^W^r) , (A9) 
L3M3 



where fc = |k| and k = k/|k|. Integrating over fi^ 3 S7^ 2 and yields 

W(i?) = 7^75 y drirf / dr 2 r|n (?7o - ri)n (r/ - r 2 ) 

x 



x 



(2tt) 3 

and the further integration over 0^ gives 



JdkiY J dk 2 Y Q( k i,k2,r ll r2,h)5i 1 i 2 5 mim2 X kl i 1 (r 1 )X k2 i 2 (r 2 ) 

h,mi i 2 ,m 2 

1 hT)2 J d 3 k(- i )^^i, 1 (A ; r 1 )i i2 (A : r 2 )jo(^)y iimi (^)y / * m2 (0 1{ ). (A10) 



W(R) = I dnr'f I di--->r 1 1 iH)(ih) - ''I )»(,(>/,. - r 2 ) 

x 



(2vr) 3 

Noting the relation (e.g., Magnus et al. 1966): 



dk t f dk 2 J2( 21 + l)Q(k u k 2 ,r 1 ,r 2 ,l)X kll (r 1 )X k2l (r 2 ) 
7 I 

' hf r f dkk'j^krMkr^MkR). (All) 



^Ur, iha-rMkR) { ^r 2 R P \ 1+ 2rlr 2 R ) (K < R < T X +r 2 ), (A12) 

(i? < |n -r 2 |, i?> n +r 2 ), 



we find 



W(-R) = 7777^-77 / / dridr 2 rir 2 n (r]Q - ri)n (ri - r 2 ) I dki I dk 2 k x k 2 



V LC 7rRJ Js 

J2(2l + l)Q(h, k 2 , n, r 2 , mhn^ik^R ~ ~ ) , (A13) 



where 5 denotes the region |n — r 2 | < R < n + r 2 . 
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So far we did not assume anything on the density field, i.e., Q(ki, k 2 , r%, r 2 , I) can be an 
arbitrary function. If we consider the case of linear theory and a scale-dependent but still local 
biasing scheme, each mode is decoupled and we can write 

Q(h, k 2 ,n,r 2 ,l) = £>i(r?o - n)£*i(% ~ r 2 )b(ki;rj - r 1 )b(k 2 ;r]o - r 2 )P(k)5(k 1 - k 2 ), (A14) 

where Di{rf) is the linear growth rate normalized to unity at present and b{k;rf) is the 
scale-dependent bias factor. In this case, equation ( |A13 ) is written as 



x J dkk 2 P(k)b(k; 7] - n)b{k;r] - r 2 ) 

„2 , Jl 7?2 



r*) 



X 

I 



£(2Z + lMknMkr^ + 2 l 2 ir2 ~ ) • (A15) 



Using the formula: 



we finally obtain 



]T(2l + lMkrMkr^ d ±A- R2 ^ = jo ( k R), (A16) 



x J dkk 2 P(k)ja(kR)b{k;r)a-n)b(k;ri -r2), (A17) 
which is identical to equation fll^). 



B. The power spectrum of mode coefficients 

Here we derive equation ( |A6| ) which is used in Appendix A. Since equation ( |A2j ) implies that 
the mode coefficient A k i m {rf) is expressed as 

Aklmiv) = J drr 2 dn^A(r],r,j)y klm (r,j), (Bl) 

one obtains 

{A klhmi (r]i)A* k2l2m2 (r] 2 )) = J t*nr?d% / dr 2 r\dVL^ 

x(A(7 ? i,ri,7i)A(r ?2 ,r 2 ,72)) i yfc 1 hm 1 (n,7i)3 ; fc 2 ; 2m2 (r2,72)- (B2) 

Since the correlation on the light-cone hypersurface does not have any special direction in 
7, the ensemble average (A(r/i, n, 7i)A(?72, r 2 , 72)) should be a function of 771, 772, J"i, ^2 and 
cos#(= 71 • 72): 

(A(r7i,n,7i)A(r/2,r 2 ,72)) = F(r] 1 ,r] 2 ,r 1 ,r 2 ,cos8). (B3) 



The function in the right-hand side of the above equation can be expanded in general as 

21 + 1 

F(rji,ri2,ri,r2,cos6) = ^ ^ C(r?i, rj 2 , n, r 2 , l)Pi(cos 6) 

i 

= E E C(r ]1 ,r]2,r l ,r2,l)Y l * m (n^ 1 )Y lm (n^). 

I m=—l 

Substituting the above expansion into equation ( [B2D and integrating over dfl^ and d£l? j2 , 
obtain equation (|A6|): 

where 

Q(h,k 2 ,m,V2,l) = / r^dri / r%dr 2 C(r)i, 772, n, r 2 , l)X kl i(ri)X k2l {r 2 ). 
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Fig. 1. — Redshift distribution of QSOs from (a) the Durham/AAT sample (Shanks & Boyle 1994) 
and that based on (b) the luminosity function by Wallington & Narayan (1993) with different 
B-band limiting magnitudes. 



-18- 
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R [h-!Mpc] 

Fig. 2. — Two-point correlation functions of QSOs for Durham/AAT sample (Shanks & Boyle 
1994; Croom & Shanks 1996) for (a) 0.3 < z < 2.2, (b) 0.3 < z < 1.4, and (c) 1.4 < z < 2.2. Our 
£,A°(R) witn M Hm = 10 11 and W 12 h~ l M Q are plotted in the shaded region. (Our Cb°( r ) is almost 
indistinguishable from rt c (i?) in this case.) For comparison, the results by Matarrese et al. (1997) 
are plotted in the dotted hatched regions. 
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Fig. 3. — Two-point correlation functions of QSOs defined on the light-cone hypersurface in cluster 
normalized standard CDM models (erg = 0.56). (a) £a(R) for the SDSS QSO catalogue for the 
threshold mass M^ m = 10 9 , 10 11 and 10 13 h~ l M & . Nonlinear mass correlation function by Peacock 
& Dodds (1996) is used for solid lines, while mass correlation function in linear theory is used 
for dashed lines; (b) the same as (a) for (c) the same as (a) for the 2dF QSO catalogue; 

(d) linear (lower curves) and nonlinear (upper curves: Peacock &; Dodds 1996) mass correlation 
functions defined on constant-time hypersurfaces z = 0, 1 and 2. 
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Fig. 4. — Dependence on the cosmological models, (a) £a(R) for the SCDM model (Qo = 1, h = 
0.5, a 8 = 0.56) of QSO for Durham/AAT sample (Shanks & Boyle 1994; Croom & Shanks 1996). 
Our £h (R) with Afi; m = 10 9 , 10 11 , and 10 13 /i _1 AfQ for bias model are plotted. Nonlinear mass 
correlation function by Peacock &: Dodds (1996) is used for solid lines, while its linear theory 
counterpart is used for dashed lines; (b) the same as (a) for £#(i?); (c) the same as (a) for LCDM 
model (r^o = 0.3, Aq = 0.7, h = 0.7, eg = 1-0); (d) the same as (c) for ^b(R)- 
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Fig. 5. — Evolution of amplitudes of two-point correlation functions at R = 15/i 1 Mpc of QSOs 
on the light-cone in (a) SCDM, and (b) LCDM models. 



